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Foreword 


As  the  accuracy  and  utility  of  numerical  ocean  models 
increases,  the  need  arises  for  realistic  atmospheric 
forcing  functions-  Since  the  wind  stress  and  heat  flux 
are  the  most  important  of  the  oceanic  forcing  func¬ 
tions  this  study  aims  to  develop  one  such  forcing 
field,  the  wind  stress,  for  the  Mediterranean  Sea.  A 
knowledge  of  the  climatological  wind  stress  fields  is 
an  essential  step  toward  determining  and  understanding 
the  wind-driven  circulation  of  any  ocean  region,  espe¬ 
cially  the  Mediterranean  Sea  where  atmosphere-ocean 
interaction  is  so  important. 


^  f)  . 

^  *  '  t* 

G.T.  Phelps,  Captain,  USN 
Commanding  Officer,  NGRDA 


Executive  Summary 


Climatological  wind  stress  averages  have  been  calcu¬ 
lated  for  the  entire  Mediterranean  Sea  for  one  degree 
latitude  by  one  degree  longitude  squares.  The  indi¬ 
vidual  stresses  were  estimated  from  ship  observations 
of  wind  speed  and  wind  direction  using  a  quadratic 
stress  law  with  variable  drag  coefficient.  Individual 
stresses  were  averaged  by  month  over  the  20-year  base 
period,  1950-1970,  to  obtain  the  climatology.  Smoothed 
versions  of  the  monthly  wind  stress  estimates,  appro¬ 
priate  for  forcing  numerical  models,  were  also 
calculated . 

Wind  stresses  exhibit  the  major  wind  patterns  that  are 
well-known  features  of  the  Mediterranean  circulation. 
The  Mistral,  Bora,  and  Etesian  wind  systems  are  the 
most  evident  of  the  area’s  features.  The  seasonal  cy¬ 
cles,  the  curl,  and  the  variance  of  the  wind  stresses 
are  presented.  Several  interesting  aspects  of  the 
relationship  of  the  Mistral  to  the  formation  of  bottom 
water  in  the  Golfe  du  Lion  region  are  discussed. 
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I.  Introduction 


Climatological  Flux  Estimates  in  the  Mediterranean  Sea 
Part  I.  Winds  and  Wind  Stresses 


The  climate  of  the  Mediterranean  Sea, 
particularly  the  wind  field,  is  perhaps 
the  most  studied  of  all  the  world’s 
oceans.  Minoan  and  Phoenecian  captains 
were  undoubtedly  aware  of  the  local  wind 
patterns  when  they  plied  the  Great  Sea 
nearly  5,000  years  ago.  This  seafaring 
tradition  was  continued,  in  turn,  with 
the  Greeks,  the  Romans,  and  the  Arabs. 

Modern  man,  while  less  dependent  on  the 
wind  as  a  means  of  propelling  commerce, 
has  found  the  interaction  of  atmosphere 
with  ocean  that  occurs  in  the  Mediter¬ 
ranean  a  compelling  subject  for  study. 
Wind-driven  surface  currents  transport 
and  diffuse  pollutants  throughout  the 
limited  volume  of  an  almost  completely 
enclosed  basin.  The  exchange  of  water 
with  the  rest  of  the  world’s  oceans  that 
does  take  place  is  largely  a  consequence 
of  the  patterns  of  evaporation  over  the 
entire  Mediterranean.  The  interactions 
of  atmosphere  and  ocean  also  result  in 
the  formation  of  deep  water  masses  in 
the  Mediterranean  that,  in  turn,  affect 
the  hydrography  of  the  world’s  oceans. 

As  part  of  the  effort  to  understand 
atmosphere-ocean  interactions  in  the 
Mediterranean  Sea,  this  exposition  fo¬ 
cuses  on  the  climatology  of  winds  and 
wind  stresses  of  the  region.  This  study 
differs  from  other  similar  works  in 
several  important  ways: 

•  Unlike  Markgraf  and  Niederschlag 
(1961)  and  comparable  weather  atlases, 
wind  stresses,  not  winds,  are  the  pri¬ 
mary  product  of  interest. 

•  This  study  uses  an  averaging  element 
that  is  1°  latitude  by  1 0  longitude,  not 
the  somewhat  standard  5°x5°  or  10°xl0°, 
or  the  irregular  regions  of  Bunker  and 
Goldsmith  (1979). 


•  These  calculations  use  more  modern  pa- 
rameterizations  than  were  available  to 
previous  authors  (e.g.,  Hellerman, 
1965). 

•  Average  stress  is  calculated  from  in¬ 
dividual  stresses,  not  from  averaged 
winds  (as  in  Hastenrath  and  Lamb,  1978). 

•  Calculations  are  performed  on  a  20- 
year  base  period,  1950-1970. 

The  primary  purpose  behind  this  effort 
was  to  develop  realistic  wind  forcing 
for  use  in  numerical  ocean  modeling. 

II.  Data  and  Methods 

The  data  used  in  the  wind  stress  calcu¬ 
lations  were  obtained  from  the  National 
Climatic  Center’s  Tape  Data  Family-11 
(TDF-11).  The  TDF-11  series  is  a  collec¬ 
tion  of  world-wide  surface  marine  obser¬ 
vations  made  by  ships  of  opportunity. 
Most  records  contain  observed  values  of 
the  quantities  necessary  for  computing 
wind  stresses  from  a  bulk  formula, 
namely:  wind  speed,  wind  direction,  air 
temperature ,  sea-surface  temperatures 
and  sea  level  pressure.  The  data  ard 
sorted  by  Marsden  squares,  eight  of 
which  cover  the  Mediterranean  Sea  (Fig. 
1).  About  1.2  million  individual  ship 
observations  spanning  the  period  1937  to 
1973  were  available  for  this  study. 

Because  of  uneven  data  coverage  of  some 
Marsden  squares,  a  20-twenty  year  period 
from  1  January  1950  to  1  January  1970 
was  selected  for  all  averaging  opera¬ 
tions.  This  period  has  relatively  uni¬ 
form  data  for  all  Marsden  squares  and 
for  all  years,  whereas  earlier  years  are 
sparsely  sampled  due  to  World  War  II, 
and  later  years  were  not  archived  com¬ 
pletely  when  the  data  was  purchased. 
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Given  the  volume  of  data  available  dur¬ 
ing  the  1950s  and  1960s,  twenty  years 
also  seems  to  be  the  minimum  time  neces¬ 
sary  to  obtain  stable  averages  on  one 
degree  grid  elements.  It  is  hoped  that 
this  will  be  a  useful  period  for  com¬ 
parison  with  other  climatological  aver¬ 
ages. 

Although  the  National  Climatic  Center 
has  edited  the  data,  further  screening 
was  performed  to  detect  and  eliminate 
erroneous  data  values.  Observations  of 
wind  speed  in  excess  of  40  meters  per 
second  (m/sec)  and  wind  directions  not 
in  the  range  from  0  to  360  were  exclu¬ 
ded  from  the  calculations.  Records  with 
missing  information  and  reported  posi¬ 
tions  not  in  or  near  the  Mediterranean 
Sea  were  ignored.  As  a  result  of  the 
editing  and  the  limited  averaging  period 
about,  800,000  ship  reports  were  used 
for  the  final  calculations. 

To  compute  wind  stresses  from  the  ship 
observations  the  well-known  nonlinear 
bulk  formula  for  turbulent  momentum 
transfer  was  used: 

•f  =  pcd  u2 

where 

t  Is  the  wind  stress, 

p  is  the  air  density, 

C^  is  the  drag  coefficient,  and 

U  is  the  wind  speed. 

In  its  most  complete  form,  p  is  a  func¬ 
tion  of  air  temperature  and  pressure, 
and  C^  is  a  function  of  wind  speed  and 
stability. 

Following  standard  practice,  the  drag 
coefficient  is  calculated  as  the  product 
of  its  value  under  neutrally  stable  con¬ 
ditions  and  a  stability  dependent  term, 

Cd(J,S)  =  F(S)  Cn(U). 

The  functional  form  of  the  neutral  drag 
coefficient  has  received  much  attention 
in  recent  years  (Garratt,  1977).  Al¬ 
though  opinion  is  far  from  unanimous. 


most  investigators  find  a  linear  rela¬ 
tionship  between  the  neutral  drag 
coefficient  and  the  wind  speed.  The 
relationship  used  for  the  wind  stress 
calculations  in  the  Mediterranean  is 

Cn(U)  =  0.9  +  0.06  U. 

This  form  was  chosen  somewhat  arbitrar¬ 
ily  as  a  compromise  among  the  several 
forms  shown  in  Figure  2.  The  scatter  of 
the  observations  is  sufficient  to  mask 
any  differences  between  the  various 
functional  forms. 

The  effects  of  stability  are  included  by 
using  a  parameterization  due  to  Rondo 
(1975)  which  uses  the  air-sea  tempera¬ 
ture  difference  and  wind  speed  to 
calculate  a  stability  parameter: 

0.1+0.03  $'+0.9  EXP  (4.8  S')  S ' <0  Stable, 

F(S  )  = 

1  .0+0,4 7  Vs*  S'>0  Unstable, 

where 

S'  =  Si  Sl/(  I  Si +0.01  Land 

S  =  (fair  -  f sea) AT  . 

Comparisons  of  stress  averages  calcula¬ 
ted  with  the  stability  correction  and 
those  calculated  with  no  correction  show 
that  stability  for  Mediterranean  Sea 
climatology  is  rather  unimportant;  dif¬ 
ferences  in  the  stresses  amounted  to 
less  than  3%. 

Although  a  variable  density,  p  ,  was  cal-( 
culated  using  the  ideal  gas  law,  compar¬ 
isons  with  constant  density  calculations 
show  that  its  effect  on  the  averaged 
stresses  was  negligible  (<1%  differen¬ 
ces)  . 

Because  of  the  nonlinear  form  of  the 
drag  law,  stress  averages  are  computed 
from  the  sum  of  individual  stresses — not 
from  averaged  winds.  Comparisons  show 
that  while  directional  differences  be¬ 
tween  the  two  techiques  are  small,  large 
(a  factor  of  2-5)  stress  magnitude  dif¬ 
ferences  do  occur.  These  magnitude  dif¬ 
ferences  are  a  reflection  of  the  wind 
variability  and  are  virtually  uncorrect- 
able . 
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Figure  2.  Comparison  of  various  functional  forms  of  the  drag  coefficient. 


The  procedure  used  to  calculate  a  vector 
average  stress  from  an  observation  of 
wind  speed  and  direction  was  to  compute 
a  stress  magnitude  from  the  wind  speed 
alone,  and  then  decompose  it  into  East- 
West  and  North-South  components  using 
the  wind  direction.  The  components  were 
then  accumulated  for  one  degree  lati¬ 
tude/longitude  "squares"  for  each  month 
of  the  year  over  the  20-year  base  per¬ 
iod.  Thus,  for  example,  the  January 
monthly  average  is  the  average  of  20 
Januarys.  The  20-year  monthly  averages 
were  used  to  form  seasonal  averages  and 
an  annual  average  as  shown  in  the  flow 
chart  of  Figure  3.  Averaging  in  this 
manner  eliminates  seasonal  biases  which 
may  develop  when  certain  seasons  have  a 
larger  number  of  data  points  than 
others . 

The  choice  of  a  spacial  averaging  inter¬ 
val  (1°  x  1°)  is  partially  based  on  a 
study  of  the  wave  number  spectrum  of  the 
wind  stress  in  the  Ionian  and  Levantine 
basins  of  the  Mediterranean.  Stresses 
were  calculated  along  ten  East-West 
lines  running  from  about  10E  to  36E, 
using  0.2°  x  0.2°  averaging  boxes.  The 
lines  were  spaced  at  0.2°  intervals  mer¬ 
idionally  starting  at  33N.  The  region  of 
coverage  is  shown  in  Figure  4.  The  ten 
series,  consisting  of  128  point  averages 
of  x  and  y  stress,  were  processed  by 
removing  the  mean  of  each  sequence, 
tapering  the  series  with  a  "ten  percent 
window,"  and  then  calculating  the  vector 
wave  number  spectrum  of  each.  The  ten 
spectra  were  then  averaged  to  give  the 
spectrum  shown  in  Figure  5. 

This  spectrum  represents  the  x-direction 
wavenumber  spectrum  of  the  averaged  wind 
stress  for  the  southeastern  Mediterran¬ 
ean  Sea.  It  shows  a  rapid  decrease  in 
energy  content  for  wavelengths  less  than 
about  250  kilometers  (km)  followed  by  a 
region  of  flattening,  which  is  interpre¬ 
ted  as  the  noise  level  of  the  averaging 
process.  While  this  wavenumber  spectrum 
is  probably  not  representative  of  the 
entire  Mediterranean,  it  does  crudely 
indicate  the  limits  of  resolution 
achievable  in  this  analysis. 


Noting  that  0.2°  boxes  are  too  small 
(energy  being  dominated  by  noise)  and 
that  5°  boxes  are  too  large  (energy 
levels  indicate  presence  of  signal)  a 
compromise  averaging  region  of  one 
degree  latitude  by  one  degree  longitude 
was  chosen  for  this  study.  This  choice 
is  further  supported  by  Saunders  (1976), 
who  found  that  wind  stress  averages  over 
regions  greater  than  that  covered  by  a 
1°  x  1°  square  resulted  in  underesti¬ 
mates  of  wind  stress  curl. 

As  with  short  temporal  averages,  small- 
scale  spacial  averages  reduce  biases 
that  can  arise  when  the  averaging  region 
has  varying  observational  densities 
coupled  with  gradients  of  the  properties 
being  averaged  (Weare  and  Strub,  1981). 
The  price  paid  for  reduced  bias  and  more 
accurate  gradient  calculation  is  lower 
stability  in  estimates  of  the  average 
due  to  the  smaller  number  of  observa¬ 
tions  for  each  averaging  interval.  In 
calculating  averages  over  regions  with 
greatly  disparate  data  densities,  the 
tradeoff  between  finer  resolution  and 
greater  statistical  stability  is  a  dif¬ 
ficult  process  to  optimize. 

Since  one  objective  of  this  work  was  to 
produce  wind  fields  suitable  for  forcing 
numerical  models,  smoothed  versions  of 
the  wind  stress  fields  have  also  been 
calculated.  Smoothing  emphasizes  the 
stability  of  the  stress  fields  rather 
than  the  resolution.  The  smoothing  oper¬ 
ation  consists  of  applying  spacial  and 
temporal  filters  to  the  monthly  averaged 
fields. 

The  temporal  filter  was  a  simple  three- 
point  ,  observation-weighted  nonrecursive 
system  of  the  form: 

F 1 (i )  «  (W  N(i-l)  +  N(i)  F(i) 

+  WN(i+l)  F(i+1))/D 

where 

W  =  0.3, 

N(i)  -  the  number  of  observations 
during  month  i 
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Figure  3.  Flow  chart  of  the  calculations 
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Figure  4.  Map  of  the  wind  stress  wavenumber  spectrum  study  region 
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Figure  5.  Vector  wavenumber  spectrum  of  the  wind  stress  with  0.2 
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F(i)  =  the  unfiltered  parameter  at 
month  i, 

FT(i)  -  the  filtered  parameter  at 
month  i, 

D  =  the  normalizing  factor,  W  N(i-l) 
+  N(i)  +  W  N(i+1). 

This  filter  reduces  the  amplitude  of  2.5 
month  fluctuations  by  half  if  the  number 
of  observations  are  the  same  during  ad¬ 
jacent  months.  If  the  number  of  observa¬ 
tions  differ,  the  filter  slightly  favors 
the  estimate  based  on  the  largest  number 
of  observations;  thus  introducing  small 
(<30°)  phase  shifts. 

The  weighting  scheme  of  the  spacial  fil¬ 
ter  is  shown  in  Figure  6.  This  filter 
also  weights  estimates  in  proportion  to 
the  number  of  observations.  It  has  cir¬ 
cular  symmetry  and  reduces  the  amplitude 
of  the  200  km  fluctuations  by  half  if 
the  number  of  observations  are  compar¬ 
able  at  all  nine  grid  elements.  This 
filter  produces  slight  spacial  phase 
shifts  when  observations  are  not  uni¬ 
form,  since  the  weights  spread  data 
information  from  well-sampled  to  poorly 
sampled  points. 

A  comparison  of  the  wind  stress  averages 
calculated  without  modification  and 
those  that  were  filtered  can  be  made  by 
comparing  Figure  7  with  Figure  12.  The 
difference  between  the  two  is  generally 
small;  however,  regions  with  strong  gra¬ 
dients  and  regions  with  limited  observa¬ 
tional  data  bases  show  significant  dif¬ 
ferences  in  the  filtered  fields.  These 
differences  are  usually  a  reduction  of 
vector  magnitude  but  occasionally  con¬ 
sist  of  vector  rotations  (<10°).  In  all 
cases  the  differences  are  not  noticeable 
over  large  regions — they  are  mainly  con¬ 
fined  to  single  points  where  the  number 
of  observations  was  significantly  smal¬ 
ler  than  adjacent  points. 

The  sources  of  error  in  climatological 
flux  estimates  have  been  discussed  by 
Bunker  (1976).  It  is  estimated  that  the 
measurements  of  wind  speed,  wind  direc¬ 
tion,  and  other  quantities  have  uncer¬ 
tainties  of  around  10%  due  to  instrumen¬ 
tal  errors  and  other  factors  such  as 


fair  weather  bias.  In  addition  the  cal¬ 
culation  of  stress  with  a  bulk  formula 
probably  introduces  uncertainties  of  not 
more  than  20%.  While  the  filtering  pro¬ 
cess  does  admit  more  uncertainty  (esti- 
imated  to  be  about  5%),  it  is  small 
compared  to  the  other  factors. 

Calculations  were  performed  on  the  Texas 
Instruments  Advanced  Scientific  Computer 
at  the  Naval  Research  Laboratory.  Care¬ 
ful  coding  and  the  vector  processing 
feature  of  this  machine  allowed  very 
efficient  handling  of  the  large  amounts 
of  data  available;  calculations  for  one 
Marsden  square  typically  used  about 
$10.00  worth  of  computer  time. 

III.  Wind  and  Wind  Stresses 

The  primary  product  of  these  calcula¬ 
tions  is  a  set  of  climatological  wind 
stress  estimates  for  each  month  of  the 
year  over  the  entire  Mediterranean  Sea. 
These  can  be  quickly  summarized  by  con¬ 
sidering  the  annual  mean  of  the  stress 
field  and  related  quantities,  and  the 
annual  cycle  as  seen  in  the  seasonal 
averages.  A  brief  description  of  these 
fields  and  the  amount  of  data  that  was 
used  is  presented  here. 

Figure  8  shows  that  the  distribution  of 
observations  roughly  follows  shipping 
lanes  in  the  Mediterranean  Sea.  The 
straits  of  Gibraltar  and  Sicily  as  well 
as  the  regions  around  the  ports  of  Mon¬ 
aco  and  Naples  have  the  highest  number 
of  observations  (10,000-12,000)  during 
the  20-year  period.  In  contrast,  the 
Adriatic  and  offshore  regions  near  Libya 
are  large  areas  where  relatively  few 
observations  (<1000)  were  made. 

The  annual  mean  wind  field  (Fig.  9)  is 
dominated  by  regional  wind  systems  hav¬ 
ing  speeds  which  range  from  5-7  m/sec. 
The  most  obvious  of  these  systems  exists 
in  the  northwestern  corner  of  the  Medit¬ 
erranean  (Golfe  du  Lion)  and  in  the 
Aegean  Sea.  Both  are  well-known  features 
of  the  region. 

Markgraf  and  Niederschlag  (1961)  have 
performed  a  similar  study  of  the  mean 
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Figure  6.  Weighting  factors  for  two  dimensional,  nine  point  spatial  filter 


10 


Q> 

U) 

m 

<D 

> 

to 

w 

w 

<y 

s- 

4-» 

W 

"D 

c 

5 


ro 

3 

C 

c 

flJ 

T3 

<u 

L- 

(1) 


LL 


V 

L. 

3 

LL 


11 


r— 

I 

o 

in 

<T> 


TJ 

O 

‘C 

<v 

o. 

<u 

£ 

4-» 

a> 

c 

‘C 

"O 

(/) 

C 

o 

-M 

tu 

> 

U 

<D 

VI 

XI 

o 


-C 

V) 


u. 

<u 

A 

E 

3 

c 


V) 

L. 

3 

O 

+-» 

c 

o 

u 


V 

i- 

□ 

a> 


12 


(/) 

"O 

c 

i 

<u 


o 

CD 

C* 

03 

i- 

<D 

> 

03 

"S 

D 

c 

c 

< 


<D 

3 

bD) 

Ll 


13 


wind  field  using  ship  observations  from 
the  period  1900-1940.  Because  their 
study  uses  comparable  data  and  the  same 
grid  (1°  x  1°),  a  comparison  can  easily 
be  made  between  their  annual  average 
(Fig.  10)  and  that  here.  The  results  are 
very  similar  despite  the  differences  in 
the  number  and  time  coverage  of  the 
observations.  Most  of  the  differences 
which  appear  in  Figure  11,  the  present 
calculations  minus  Markgraf’s  data,  are 
quite  small  with  respect  to  the  mean 
wind  field.  The  larger  deviations  can 
easily  be  explained  as  a  result  of  sta¬ 
tistical  uncertainty  in  the  averages. 
Systematic  differences,  which  can  be 
seen  along  the  southern  coastlines  of 
the  Mediterranean,  are  probably  too 
small  to  attribute  to  long-term  secular 
change  and  may  be  the  result  of  some  of 
the  errors  discussed  above. 

The  annual  mean  wind  stress  (Fig.  12) 
naturally  exhibits  patterns  similar  to 
that  of  the  wind  field.  Its  major  fea¬ 
ture  is  the  strong  southeastward  stress 
in  the  western  Mediterranean.  Uniform 
eastward  stresses  in  the  Tyrhennian  Sea 
contrast  with  westward  stresses  in  the 
northern  reaches  of  the  Adriatic.  Stress 
patterns  in  the  eastern  Mediterranean 
largely  consist  of  eastward  to  south¬ 
eastward  stresses  in  the  Ionian  and  Lev¬ 
antine  basins,  and  southward  stresses  in 
the  Aegean  Sea.  Magnitudes  range  from 
highs  of  nearly  3  dynes/cm^  (0.3  pas¬ 
cal)  in  the  northwestern  Mediterranean 
to  0. 5-1.0  dyne/cm^  over  much  of  the 
eastern  basin. 

The  annual  mean  wind  stress  magnitude, 
contoured  in  Figure  13,  has  values  that 
range  from  almost  3.0  dynes/cm 2  in  the 
western  Mediterranean  to  less  than  0.8 
dyne/cm^  in  large  regions  of  the 
eastern  basin  and  in  the  Alboran 
Sea/Baleric  Sea  region  of  the  western 
basin.  The  Aegean  Sea  appears  to  be 
another  area  of  relatively  strong  winds 
with  magnitudes  around  1.2  dynes/cm^. 

The  vector  curl  of  the  wind  stress  is  an 
important  quantity  which  is  readily 
calculated  from  the  wind  stress  field. 


Figure  14  shows  the  curl  calculated  from 
the  annual  mean  wind  stress  using  cen¬ 
tered  differences  to  estimate  deriva¬ 
tives.  The  curl’s  importance  stems  from 
its  relation  to  the  divergence  of  wind- 
driven  surface  currents  and,  hence,  from 
a  proportionality  to  vertical  velocity 
at  the  base  of  the  so-called  Ekman  layer 
(see  Stommel,  1965).  Given  the  short 
north-south  extent  of  the  Mediterranean 
and  the  relatively  minor  variations  of 
the  Coriolis  parameter  that  this  en¬ 
tails,  the  curl  fields  should  be  good 
indicators  of  vertical  velocity  with  the 
rough  conversion  that  curl  of  1x10”® 
dynes/cm^  is  equivalent  to  a  vertical 
velocity  of  10  cm/day.  Two  regions  of 
significant  curl  appear  in  the  figure: 
the  Western  Basin  and  the  Aegean  Sea. 
The  magnitude  of  the  mean  annual  verti¬ 
cal  velocities  that  can  be  inferred  from 
the  stress  curl  to  south  of  France  is 
from  10  to  50  cm/day:  upwelling  in  re¬ 
gions  of  positive  curl  and  downwelling 
in  regions  of  negative  curl.  Vertical 
velocities  in  the  eastern  Aegean  are 
almost  comparable . 

Contours  of  mean  annual  wind  stress  var¬ 
iance  (Fig.  15)  show  that  wind  stress  in 
the  western  Mediterranean  is  highly  var¬ 
iable  as  well  as  strong.  The  rest  of  the 
Mediterranean  exhibits  relatively  uni¬ 
form  low  variability  even  in  regions 
where  the  winds  are  strong,  such  as  the 
Aegean  Sea. 

On  a  seasonal  basis,  spring  (March-May) 
stresses  are  characterized  by  moderate 
(2  dynes/cm^)  southwestward  wind 
stress  in  the  western  Mediterranean  and 
weak  stress  in  the  Aegean  Sea  and  Ionian 
basin.  The  mean  spring  stress  field, 
shown  in  Figure  16,  also  shows  the  east¬ 
ward  to  southeastward  stress  of  about 
0.5  dyne/cm^  that  characterizes  the 
Levantine  basin  throughout  much  of  the 
year.  Another  point  to  note  is  the  pres¬ 
ence  of  weak  southeastward  stress  along 
the  Algerian  coast. 

Summertime  stresses  (Fig.  17)  show 
weakening  (1  dyne/cm^)  in  the  western 
Mediterranean  but  strengthening  in  the 
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Figure  14.  Contoured  annual  wind  stress  curl 
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Aegean  Sea  and  Levantine  basin.  Stresses 
along  the  Algerian  coast  are  southward 
to  southwestward  in  contrast  to  the 
springtime  stresses . 

Autumn  stress  fields  (Fig.  18)  are  very 
similar  in  strength  and  pattern  to  those 
of  summer.  Exceptions  occur  in  the  Ion¬ 
ian  basin  where  southward  stress  shifts 
to  weak  eastward  stress,  and  along  the 
Algerian  coast  where  stress  reverses 
again  to  be  come  eastward. 

Not  surprisingly,  winter  stresses  (Fig. 
19)  are  generally  strong  over  the  entire 
Mediterranean,  with  the  exception  of  the 
Aegean  Sea  where  the  stress  is  weak  and 
appears  to  be  controlled  by  a  small  cy¬ 
clonic  feature.  The  entire  eastern  Med¬ 
iterranean  is  dominated  by  relatively 
strong  (1  dyne/cm^)  eastward  stress 
while  the  northwestern  Mediterranean  is 
dominated  by  2—3  dynes/cm3  southeast¬ 
ward  stress.  Several  other  regions  also 
show  significant  increases  in  averaged 
stress.  The  northern  reaches  of  the 
Adriatic  Sea  have  strengthened  stresses 
to  the  southwest  and  the  Alboran  Sea 
shows  freshened  stresses  to  the  east. 
Because  of  its  strength,  the  winter 
stresses  are  clearly  the  most  important 
contribution  to  the  annual  average. 

Monthly  mean  stresses  are  presented  in 
an  appendix.  They  exhibit  the  seasonal 
tendencies  described  above  and  show  the 
transitions  from  season  to  season  that 
justify  the  divisions.  From  these  month¬ 
ly  means  it  is  possible  to  assemble  a 
time  series  of  vector  stresses  for  a 
particular  region  by  averaging  the 
stresses  over  the  area  for  each  month 
and  plotting  the  twelve  averages  on  a 
time  axis.  Figures  20-25  show  the  yearly 
cycle  of  wind  stress  for  various  Medi¬ 
terranean  basins. 

These  figures  reiterate  most  of  the 
points  inferred  from  the  seasonal  stress 
plots.  As  noted  above,  the  northwestern 
Mediterranean  has  the  greatest  stresses, 
with  a  maximum  magnitude  in  December  and 
a  minimum  in  September.  The  direction  of 
the  stress  is  consistently  around  145°. 


By  contrast,  the  Aegean  Sea  stress  ser¬ 
ies  shows  a  maximum  in  late  summer  or 
early  autumn  and  a  minimum  during  the 
winter  months.  The  stresses  in  the  Lev¬ 
antine  basin  are  the  most  consistent 
throughout  the  year  with  magnitudes  be¬ 
tween  0.3  and  0.6  dyne/cm^  and  direc¬ 
tions  that  range  from  100°  to  125°.  The 
series  in  the  southwestern  Mediterranean 
(south  of  the  Baleric  Islands)  clearly 
shows  the  stress  reversal  mentioned 
above « 

An  interesting  additional  calculation, 
which  can  be  performed  with  the  annual 
averages  of  winds  and  wind  stresses,  is 
the  regression  of  the  averaged  stress  on 
the  averaged  wind  squared.  Figures  26 
and  27  show  the  scatter  plots  of  these 
quantities  for  both  the  western  and 
eastern  portions  of  the  Mediterranean 
and  the  least-squares  regression  lines. 
The  slope  of  the  regression  lines  can  be 
related  to  a  drag  coefficient  for  use 
with  averaged  winds.  These  slopes  indi¬ 
cate  a  drag  coefficient  of  5.6  x  1CT3 
for  the  western  basin  and  2.3  x  10~3 
for  the  eastern  basin.  The  difference 
between  the  two  drag  coefficients  is 
indicative  of  the  difference  in  wind 
variability  for  the  eastern  and  western 
basins  of  the  Mediterranean.  The  differ¬ 
ence  also  justifies  the  calculation  and 
averaging  of  individual  stresses,  since 
a  priori  knowledge  of  drag  coefficients 
for  averaged  winds  was  not  available. 
The  conclusion  is  that  if  the  averaged 
winds  had  been  used  with  a  typical  drag 
coefficient  of  1.5  x  10“3,  the  stres¬ 
ses  would  have  been  severely  underesti¬ 
mated  . 

IV.  Discussion 

A  number  of  the  wind  systems  that  cause 
the  wind  stress  patterns  noted  are  so 
persistent  that  they  have  been  given 
names:  Mistral,  Bora,  Scirocco,  Levant¬ 
er,  Vendaval,  Libeccio,  and  Etesians  are 
names  of  several  regional  wind  systems 
in  the  Mediterranean.  Only  the  strongest 
of  these  systems  (the  Mistral,  the 
Etesians,  and  the  Bora)  were  detectable 
in  the  climatological  averages,  however. 
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ure  19.  Winter  wind  stress  average  for  points  with  more 
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Figure  20.  Year-long  cycle  of  stresses  in  the  northwestern  Mediterranean 


SOUTAHESTERN  MEDITERRANEAN 


(2**W3/S3NXQ)  5S3yiS 


27 


Figure  21.  Year-long  cycle  of  stresses  in  the  southwestern  Mediterranean 
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Figure  22.  Year-long  cycle  of  stresses  in  the  Tyrhennian  Sea 
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Figure  23.  Year-long  cycle  of  stresses  in  the  Ionian  Sea 
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Figure  24.  Year-long  cycle  of  stresses  in  the  Aegean  Sea 
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Figure  25.  Year-long  cycle  of  stresses  in  the  eastern  Mediterranean 
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western  Mediterranean 


A.  The  Mistral 

The  Mistral  is  the  intense  northwesterly- 
wind  that  appears  during  all  seasons  in 
the  western  basin  of  the  Mediterranean. 
It  blows  down  Rhone  River  valley  through 
the  gap  between  the  Pyrenees  and  the 
Alps  (the  Garronne-Carcassone  gap).  It 
is  a  cold,  dry,  and  very  persistent  wind 
that  meets  the  Mediterranean  at  the 
Golfe  du  Lion  and  is  responsible  for  the 
large  southeastward  stresses  in  the  cli¬ 
matological  averages.  The  heat  loss 
which  results  as  the  Mistral  sweeps 
across  the  northwestern  Mediterranean 
during  the  winter  is  thought  to  be 
responsible  for  the  large  regions  of 
convective  overturning  that  have  been 
observed  in  the  area  (MEDOC  Group, 
1970).  The  Mistral  is  also  responsible 
for  the  large  values  of  wind-stress  curl 
that  have  been  calculated  from  the 
stress  averages.  Its  narrow,  jet-like 
profile  and  strong,  persistent  stress 
magnitudes  result  in  large  gradients 
across  the  axis  of  the  wind  and,  hence, 
in  large  wind-stress  curls;  negative 
curl  to  the  right  of  the  Mistral’s  axis 
and  positive  curl  to  the  left  of  the 
axis  looking  downwind. 

The  vertical  motion  that  results  from 
the  divergence  of  horizontal  currents  in 
the  Ekman  layer  Is  substantial  for 
oceanic  conditions.  Vertical  velocities 
of  up  to  70  meters/day  are  possible, 
given  the  magnitude  of  the  curl  computed 
in  the  Golfe  du  Lion  region.  Since  these 
figures  are  long-term  averages,  individ¬ 
ual  months  or  seasons  are  liable  to 
exhibit  even  larger  vertical  movements. 
The  effects  on  the  mean  circulation  in 
the  region  should  certainly  be 
substantial  enough  to  observe  with 
classical  oceanographic  surveys . 

Hydrographic  surveys  of  the  Northwestern 
Mediterranean  do  show  the  presence  of  a 
large  cyclonic  gyre  centered  approxi¬ 
mately  at  42N  and  5E  (Ovchinnikov, 
1966).  This  circulation  feature  is  a 
reasonable  consequence  of  the  positive 
curl  and  upward  motion  that  is  inferred 
from  the  wind  stress  calculations. 


Direct  measurements  of  vertical  velocity 
by  Voorhis  and  Webb  (1970)  during  the 
MEDOC  experiment  indicate  that  large  (2 
cm/sec)  upward  velocities  do  exist  in 
the  region  and  are  persistent  over  per¬ 
iods  of  several  days.  Calculations  of 
vertical  velocities  by  Seung  (1980)  also 
show  substantial  upward  movement  on 
monthly  time  scales.  More  significantly, 
Swallow  and  Canton  (1973)  suggest  that 
the  gyre  plays  an  important  role  in  the 
initial  or  preconditioning  stages  of 
bottom  water  formation  by  bringing  the 
salty  Intermediate  water  close  to  the 
surface.  The  possibility  that  the  Mis¬ 
tral  Is  responsible  for  preparing  the 
northwestern  Mediterranean  for  bottom 
water  formation  as  well  as  for  the 
actual  cooling  of  the  surface  water  is 
an  interesting  topic  for  future  study. 

Corroborating  evidence  for  the  large 
downward  velocities,  which  should  occur 
to  the  west  of  the  Mistral’s  axis,  has 
not  been  found.  Hydrographic  surveys  do 
not  seem  to  show  any  indication  of  an 
anticyclonic  gyre.  Measurements  that 
indicate  downward  bursts  of  water  during 
the  intense  cooling  of  winter  are  clear¬ 
ly  irrelevant  to  the  Ekman  layer  argu¬ 
ments  advanced  here,  but  it  should  be 
noted  that  the  direct  measurements  cited 
are  all  to  the  east  of  the  region  of 
negative  wind-stress  curl. 

B.  The  Etesians 

In  the  eastern  Mediterranean  the  dom¬ 
inant  winds  are  the  Etesians  (Greek)  or 
Meltems  (Turkish).  These  winds  are  fun- 
neled  on  to  the  Aegean  Sea  through  the 
gap  that  exists  between  the  mountains  of 
the  Balkans  and  the  mountains  of  Ana¬ 
tolia — i.e.,  along  the  Rosporos  and 
Dardanelles,  and  the  Maramara  Sea.  The 
Etesians  are  also  northerly  winds,  but 
are  strongest  in  the  late  summer  or 
early  autumn.  The  region  of  positive 
wind-stress  curl  which  hugs  the  Turkish 
coast  Is  another  area  that  has  been 
indicated  in  the  formation  of  subsurface 
water  masses  in  the  Mediterranean  (Wust, 
1961).  There  seems  to  be  no  evidence, 
however,  that  the  Etesian  wind  system 
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has  any  effect  on  the  hydrography  of  the 
area  or  the  process  of  cooling. 

C  The  Bora 

The  last  of  the  Mediterranean  wind  sys¬ 
tems  to  be  evident  in  the  20-year  aver¬ 
ages  is  the  Bora.  The  Bora  is  another 
wind  that  is  directed  by  the  topography 
of  the  region — in  this  case,  the  Trieste 
gap  at  the  northern  end  of  the  Adriatic 
Sea.  This  break  between  the  European  and 
Dinaric  Alps  directs  a  westward  flow 
over  the  northern  Adriatic.  The  Bora  is 
a  strong  but  infrequent  winter  wind; 
consequently,  it  appears  only  as  a 
slightly  larger  westward  stress  in  the 
winter  averages . 

V.  Summary 

Twenty-year  mean  winds  and  wind  stresses 
have  been  calculated  from  the  ship 
observation  reports  of  the  National  Cli¬ 
matic  Center.  The  winds,  stresses,  and 
derived  quantitites  are  organized  into 
monthly,  seasonal,  and  annual  mean 
fields  with  one  degree  resolution. 
Smoothed  versions  of  all  fields  were 
also  prepared.  Estimates  of  the  wind 
stress  were  obtained  using  an  aerodynam¬ 
ic  drag  law  to  calculated  stresses  from 
individual  ship  reports.  The  stress 
estimates  were  then  averaged  on  the  one 
degree  squares  for  each  month  for  the 
20-year  period. 

The  average  wind  and  stress  field  show 
the  major  wind  systems  of  the  Mediter¬ 
ranean  Sea.  Of  these  systems,  the  Mis¬ 
tral  and  Etesian  winds  are  the  most 
obvious.  The  Mistral  is  of  particular 
interest  due  to  its  influence  on  bottom 
water  formation  and  wind  driven  circula¬ 
tion  in  the  western  Mediterranean. 
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Appendix:  Monthly  Mean  Stresses 
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Figure  A2.  February  wind  stress  (filtered) 
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Figure  A3.  March  wind  stress  (filtered) 
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